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Abstract 

In universes with significant curvature or cosmological constant, cosmic microwave 
background (CMB) anisotropies are created very recently via the Rees-Sciama or in- 
tegrated Sachs- Wolfe effects. This causes the CMB anisotropies to become partially 
correlated with the local matter density (z < 4). We examine the prospects of using 
the hard (2-10 keV) X-ray background as a probe of the local density and the mea- 
sured correlation between the HEAOl A2 X-ray survey and the 4- year COBE-DMR 
map to obtain a constraint on the cosmological constant. The 95% confidence level 
upper limit on the cosmological constant is ^ 0.5, assuming that the observed 
fluctuations in the X-ray map result entirely from large scale structure. (This would 
also imply that the X-rays trace matter with a bias factor of 6a; ~ 5.6 .) This 
bound is weakened considerably if a large portion of the X-ray fluctuations arise 
from Poisson noise from unresolved sources. For example, if one assumes that the 
X-ray bias is bx = 2., then the 95% confidence level upper limit is weaker, J7a ^ 0.7. 
More stringent limits should be attainable with data from the next generation of 
CMB and X-ray background maps. 
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The recent observations of anisotropics in the microwave background (Bennett, 
et al., 1996a; Netterfield et al., 1997; Gundcrson, et al., 1995; Scott, Silk & 
White, 1995) have been widely interpreted as the result of potential, tem- 
perature or velocity fluctuations at the surface of last scattering, originating 
at very high red shifts {z > 1000). However, in many cosmological models, 
a significant portion of the anisotropies are produced much more recently. In 
critical density cosmologies, nonlinear gravitational clustering produces a time 
dependent gravitational potential, as bound objects turn around and collapse. 
This gives rise to CMB anisotropies on small angular scales. In cosmologies 
with less than critical density (e.g. flat universes with a large cosmological 
constant and open universes), time dependence of the gravitational potential 
is induced by a change in the expansion law of the universe at late times, 
even in the linear perturbation regime. This gives rise to CMB fluctuations 
on larger scales. The nonlinear effect is usually referred to as the Rees-Sciama 
(RS) effect (Rees & Sciama, 1968), and the linear effect as the integrated 
Sachs- Wolfe (ISW) effect (Sachs & Wolfe, 1967). An earlier paper by two of 
us (Crittenden & Turok, 1996) referred to both as the Rees-Sciama effect. 

The recently produced CMB fluctuations result from time variations in the 
gravitational potential and are correlated with the nearby matter density {z < 
4). Observing such correlations could place a strong constraint on important 
cosmological parameters, such as the matter density, Q„i, and the cosmological 
constant, Qa- Above all, such correlations would offer a rare opportunity to 
observe CMB anisotropies as they are being produced (Crittenden & Turok, 
1996). 

Implementing this, however, requires a probe of the matter density at high 
red shifts. Possible probes include radio galaxies and quasars, and a number 
of large scale surveys of these objects are currently underway. In this paper, 
a map of the hard X-ray background is used as a measure of the local matter 
density. X-rays with energies of order a few keV appear to be produced pri- 
marily by active galactic nuclei (AGN) (Comastri et al., 1995) and so should 
reflect the mass distribution on large scales. Here, we investigate the cosmo- 
logical limits which result from cross correlating the HEAOl A2 2-10 keV 
X-ray map (Boldt, 1987) with the four year GOBE DMR map of the cosmic 
microwave background (Bennett, et al., 1996a). 



1 X-ray emission as a tracer of mass 

In order to use the X-ray background in the type of analysis suggested above 
two questions must be answered: "How well do X-ray sources trace mat- 
ter?" and 'What is the redshift distribution of the X-ray emission?". The 
latter of these questions will be addressed using a particular model of X-ray 



2 



sources (Comastri et al., 1995) to compute the intensity-rcdshift distribution, 
di{z)/dz, i.e. the portion of the intensity of the X-ray background that arises 
at redshift z. The resuhs presented in this paper are not very sensitive to 
the particular form of di{z)/dz as long as a signficant portion of the intensity 
arises from redshifts around ~ 1. This is the case for most current models of 
the X-ray background. How well X-rays trace matter is less well known. In this 
paper we will assume a simple linear bias, i.e. that the fractional fluctuations 
in X-ray emissivity are linearly proportional to the fractional mass density 
fluctuations, {5ex/ex) = bx{Spm/ Pm), where is the bias factor. The linear 
bias assumption would hold if the X-ray sources were produced at the high 
peaks of a Gaussian random field, in which case the bias would be approx- 
imately the height u of the peaks, in units of the standard deviation on the 
appropriate smoothing scale. More generally, it is quite plausible that X-ray 
sources trace the matter density, and the linear bias model may be taken as a 
crude parametrization of this. For example, if X-ray sources are produced by 
local physical processes which are completely independent of the long wave- 
length density perturbations, one would have bx — 1, whatever the statistics 
of the density field. 

While the source of the X-ray background is still very much an area of intense 
research, it is already clear that discrete sources are the major component. 
Deep ROSAT observations have resolved > 60% of the soft (0.5 - 2.0 keV) X- 
ray background into discrete sources (Hasinger et al., 1993) and unified AGN 
(Seyfert galaxies and quasars) models have successfully described both the 
spectrum and the brightness of the hard (2-10 keV) background (Comastri et 
al., 1995; Madau, et al., 1994; Zadarski, Zycki & Krolik, 1993) as well as the 
local volume emissivity (Miyaji, et al., 1994). Despite the successes of these 
models, it is quite possible that they will have to be modified in some way 
(e.g. more detailed evolution, refined spectra, a less 'unified' AGN model, a 
new population of sources, etc.) as more is learned about X-ray sources. In 
particular, correlation and fluctuation analyses are still not well understood in 
either the hard or soft bands (Miyaji, et al., 1994; Chen et al., 1994; Carrera 
et al., 1993). However, it seems unlikely that the average flux distribution, 
di{z)/dz, will be greatly changed by such modiflcations. Given that the ISW 
effect is not particularly sensitive to the detailed shape of di{z)/dz, we do not 
believe that such modifications to the model will have a large effect on the 
results of this paper. 

The estimate of the X-ray bias factor, b^, is more problematic. Measurements 
of the local X-ray bias are still uncertain. Prom the dipole anisotropy of the 
local AGN distribution, Miyaji (1994) estimated that b^^;;;^^'^ ~ 3.5// where 
/ is fraction of the total gravitational acceleration of the Local Group con- 
tributed by matter within 45 h^^ Mpc, Qm is the matter density, and h is the 
Hubble constant in units of 100 km s~^ Mpc~^. For flrn = 0.3 and / ~ 0.5 this 
implies b^ — 3. However, it is not clear if the bias measured in such nearby 



3 



X-ray sources is appropriate on the much larger scales of interest for the ISW 
effect. For example, it is possible that emission from nearby clusters might 
make the local bias appear higher. 

In the context of a given theory, where the power spectrum of density fluc- 
tuations is known and appropriately normalized by the COBE detection of 
the anisotropics, it is also possible to place an upper limit on the X-ray bias 
by observing the angular auto-correlation function (ACF) of the X-ray back- 
ground. One can compare this to that predicted by the theory assuming a 
linear, scale- independent bias. Unfortunately, the large scale (> 1°) angular 
correlations of the X-ray background have not been accurately measured. From 
our own measurements of the ACF of the hard background we found that, in 
the context of cold dark matter (CDM) models with a cosmological constant, 
bx ^ 5.6f^^^^ (see Table 1). This should be considered as an upper limit since 
the HEAO beam size is large (3° x 3°) and can itself induce correlations even 
when viewing randomly placed sources. While limits found in this manner are 
necessarily model dependent, they do provide important context for assessing 
the viability of a model. 

One might argue by analogy with other populations. Ordinary galaxies have 
correlation lengths of order Tq — 5h^^ Mpc (Groth & Peebles, 1977) and 
are thought to be only mildly biased, 6^ ~ 1. This might be considered a 
reasonable lower limit to the X-ray bias. The correlation lengths, of radio 
galaxies (Loan, Wall, & Lahav, 1997; Sicotte, 1995), QSO's (Andreani & Cris- 
tiani, 1992), and groups of galaxies (Bahcall, 1996) are all on the order of 10 
- 18h~^ Mpc. If, as seems likely, hard X-ray sources have similar clustering 
properties to these, then we expect the X-ray bias to be of order — 2 — 3. 
On the other hand, rich clusters of galaxies are even more strongly correlated 
(Bahcall, 1996), Vg — 20 — 25/i^^Mpc and if X-ray sources are instead like these 
it would imply b^ > 3. It is well known that the X-ray background is highly cor- 
related with these other objects (Rcfregier, Helfand & McMahon, 1997; Soltan 
et al., 1996; Soltan et al., 1997; Miyaji, et al., 1994); however, the correlation 
length is still uncertain. It is clear that the X-ray bias (if indeed linear bias 
is a good description of the distribution of X-ray sources) is not well known. 
Therefore, we have chosen to report our results with X-ray bias as an unknown 
parameter. 

Finally, as an independent, qualitative check on how X-rays trace mass we have 
cross-correlated the HEAO 2-10 keV map with the appropriately smoothed 
and binned Greenbank 5 GHz radio source counts (Gregory & Condon, 1991). 
The correlation coefficient of these two maps is signficant, {SnSi) / Sirms^iT'rms — 
0.2 dz 0.05, indicating that they are tracing correlated populations. The cross- 
correlation function can also be expressed as {dndi) ~ {^r) + WRx{n){Lx) 
where lr is the X-ray intensity of the radio sources, ix the intensity of other 
X-ray sources in the fleld, {n) is the average number of radio sources in a 



4 



resolution element, and wxr is the cross-correlation coefficient between radio 
sources and the X-ray sources. Since the numerical value of {dndi) / (l) is 0.03, 
either 3% of the X-ray background is due to the 5 GHz radio sources or the 
correlation coefficient between radio sources and X-ray sources is on the or- 
der of Wxr — 10~^. In either case, the imphcation is that X-ray sources are 
signficantly clustered. The Greenbank radio sources are themselves strongly 
correlated with correlation lengths of 10 — 18h^^ Mpc (Loan, Wall, & La- 
hav, 1997; Sicotte, 1995). When smoothed on the HEAO beam scale the radio 
source correlation coefficient is wrr ~ 2 x 10~^. These values are comparable 
to wxx deduced from the X-ray auto-correlation function (see Section 3.1). 
The Greenbank survey is dominated by moderate redshift (^ ~ 1) sources 
(Loan, Wall, & Lahav, 1997; Sicotte, 1995) which implies (with some assump- 
tion about X-ray source evolution) that most of contribution to the cross- 
correlation, {5n5i), comes from the redshift range 0.1 < 2; < 1. This redshift 
interval is the source of rougly 25% of the X-ray background (see Fig. 1). Simi- 
lar results were found cross-correlating the X-ray background with the FIRST 
radio survey (Cress et al., 1996). These results cannot be made more quanti- 
tative without a detailed analysis of the X-ray/radio source cross-correlation 
which is beyond the scope of this paper. 



2 CMB and X-ray Correlations 

We will focus here on the two point correlations of the microwave temperature 
and X-ray intensity, and the related power spectrum of the density fluctua- 
tions. The density can be written as 



where p{t) is the average matter density at proper time t. The power spectrum 
is defined as the Fourier transform of the two point correlation function and 
is given by, (f^k^k/) = Pfe5(k — k') where 5(x, to) = Z^k^kC**^'^ and to is the 
present time. (Note that we will use x and k to refer to comoving position 
and wavenumber.) 

All sky maps of the X-ray intensity or temperature fluctuations can be nat- 
urally expanded in terms of spherical harmonics. For example, we can write 
the fluctuation in X-ray intensity as 



p(x,t) =p(t)(l + 5(x,t)). 



(1) 



(2) 
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where 5/,(n) = /,(n) — I and n is the unit direction vector. A similar expression 
can be written for the temperature anisotropy 5T(n). The expectation value 
of the moments defines an angular power spectrum 

Q = (|ai^|2), (3) 
which is related to the expectation value of the autocorrelation function by 
'5t(n) 5i{-n')\ 



^^Ei^l + l)CiPi{cose), (4) 



47r ^ 



where cos^ = n • n' and Pi{x) is a Legendre polynominal. 

2.1 CMB Anisotropies 



In the approximation of instantaneous recombination, the microwave anisotropy 
in a direction n on the sky is given in Newtonian gauge by 

— (n) = [-S, + V • n + + 2 1 dr^r, n(ro - r)). (5) 

Tr 



The integral is over the conformal time r, and Tr and Tq are the times of 
recombination and the present, respectively. The first term {\S^) represents 
density perturbations of the radiation-baryon fiuid, the second term is the 
Doppler shift (v • n), and the third is the Newtonian potential ($), where all 
of these are evahiated on the surface of last scattering. The last term, usually 
called the Integrated Sachs- Wolfe (ISW) term, represents the effect of a time 
varying gravitional potential along the line of sight. Heuristically, it represents 
the redshifting of photons which 'climb out' of a different potential than they 
'fell into'. 

The CMB anisotropies are the sum of contributions created near the surface 
of last scattering and those created recently by the ISW effect. Since the two 
contributions are associated with perturbations of very different wavelengths, 
they are nearly uncorrelated except on the largest angular scales. The moments 
can be written as. 
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and the fact that they are weakly correlated implies 



(7) 



where C/^^ = (|a[^^n and Ct' = (lafip). 

The ISW term arises when the CMB photons pass through a potential fluc- 
tuation which is evolving in time, i.e. 

af^^ = 4m' E ^ / 2^(^, r)ji{kr)drYi^{n^), (8) 

k -ffe 

where r is the conformal distance, tq — r, and ji (x) is a spherical Bessel func- 
tion. Squaring this and integrating over the directions of k, we find 

(9) 

2.2 The X-ray Model 

To evaluate the expected X-ray brightness fluctuations, a model is needed 
which describes the redshift distribution of the sources, as well as their lumi- 
nousities. We consider the unified AGN model of Comastri et al. (1995), which 
reproduces source number counts and the flux and spectrum of the X-ray back- 
ground in both soft and hard X-ray bands. In this model, AGN are divided 
into two luminosity classes, Scyfcrt galaxies and quasars, which range in X- 
ray luminosity from 10''^ ergs s^^ to 10"^^ ergs s"-^ (0.3-3.5 keV). The power-law 
luminosity function has a break at a luminosity oi Lb = 6.9 x 10"^^ ergs s~^ 
where AGN with luminosities L < Lb are designated Seyferts and those with 
L > Lb are designated quasars. 

The AGN are thought to be surrounded by an absorbing molecular torus, with 
the amount of absorption depending on the line of sight. As an approximation, 
the AGN are divided into five absorption classes characterized by their H col- 
umn density, with column densities ranging from zero up to Nh = 10^''cm~^. 
The number density of each class is proportional to the number density of 
unabsorbed AGN, i.e. ns,i = Wins,o and j = wiUq^q where nsfl (^q,o) is the 
total number of unabsorbed Seyferts (quasars) and the index i runs over ab- 
sorption class. The relative density in each class is assumed to be independent 
of both source luminosity and redshift. A detailed description of this division 
and the associated spectra can be found in Comastri et al. (1995). 

The average number of AGN in a volume at redshift z is dN{z) = n{z)dV 



CI 



ISW 



kMk 



2^{k,r)ji{kr)dr 
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where dV = d\{z)dD,dl, dA^z) is the angular diameter distance, dD, is the sohd 
angle and dl is proper radial distance. Proper distance can be expressed as 

, da dz a 

dl = dt = — = , 10 

a' 1 + za' ^ ^ 



where "a" is the expansion factor and primes indicate derivatives with re- 
spect to proper time. In the Comastri model, the comoving number density 
is assumed to be constant, so that n{z) = nQ{l + z)^ where Uq is the to- 
tal unabsorbed AGN density. It can be expressed as, Hq = us^ + uq^q = 
■^imin 'Ps{L)dL + /^^"'^ (t)q{L)dL, where 0s is the luminosity function of Seyferts 
and 0q the luminosity function of quasars. 

Another key feature of most X-ray models is that the X-ray luminosity evolves 
with redshift as 

L(^) =L(0) X {l + zf. (11) 



The Comastri model assumes (5 — 2.Q and a cutoff redshift, z cutoff — 2.25, 
above which the luminosity is assumed to be constant, up to maximum redshift 
of Zmax = 4. There is some evidence the evolution could be even stronger 
than assumed by Comastri et al., with a slightly higher index, /5 — 3, and a 
somewhat lower cutoff redshift, Zcutoff — 1-8. A detailed examination of these 
luminosity models is found in Boyle et al. (1994). 

The flux received from an object at a given redshift is related to its luminosity 

by 



where the luminosity distance rf^ is related to the angular distance by dL — 
dA{l + zY . For band limited detectors this relation becomes 



f{z, I)L{z) 



^(.,L,/)^^^^. (13) 



where L{z) now represents the energy emitted in the detector band [Ei,E2\ 
and f{z,I) is the usual K-correction 



£/ T\ ^detected _ !E:i^:i mdE 
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_ (1 + IjEjl + z))dE 

Since the K-correction depends on I{E), the spectrum of the source, it is a 
different function for each absorbtion class. 

At a given redshift the flux averaged over the various types of AGN is obtained 
by summing over both luminosity and absorption, i.e.. 



^^^'^^ ^{Ls)j:^if{z,Isi) + {LQ)Y.^J{z,lQ,)] (15) 



indl 



where (Ls) = ji:^^^(t>s{L)LdL/noY.Wi and {Lq) = J^-^ (f)g{L)LdL/noEwi. 
We denote the term in parentheses as Lk{z), the average K-corrected lumi- 
nosity per unit volume. For the absorbed spectra used by Comastri et al., this 
factor can be approximated by Lk{z) oc (1 + ^)/(1 + 0.15^;). 

Finally, the "space weighting function", i.e. the contribution to the band lim- 
ited X-ray intensity as a function of redshift, is 

^?M^W^„.(i + .)=dida^^(i±fL.(.). (16) 



Using the fact that a'/a — -f^o(^A + ^m(l + z)^y^'^ for the flat cosmologies 
studied here, and simplifying, the weighting function is 

dNjz) ^ no{l + z)^-''LK{z)dQ 

dz A7rHoinA + n^{i + z)^y/^' ^ ' 



Figure 1 is a plot of dN{z)/dzT{z) in arbitrary units. It is clear from this 
figure that much of the flux arises from sources at redshifts of order z ~ 
1 — 2 as required in order to detect a significant ISW effect. The shape of 
the distribution is fairly flat and is similar to (though it extends to somewhat 
higher redshifts) the distribution used in Crittenden & Turok (1996) which 
used dN{z)/dz measured in a flux limited survey and assumed that the flux 
received from an object at a given redshift was constant. 

It should be noted that for their analysis Comastri et al. (1995) adopted a 
Hubble constant of Hq = 50kms^^Mpc^^ and a deceleration parameter of 
qo = 0. The fact that these differ from the cosmological parameters of the 
models in this paper is not of great concern since ultimately we only use the 
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Fig. 1. The redshift distribution of contributions to the intensity fluctuations in the 
cosmic X-ray background, based on the model of Comastri, et al. The cutoff at low 

redshift results from excluding high flux X-ray sources, while the cusp at z = 2.25 
results from the change in the luminosity evolution of the sources. 

space weighting function which is directly tied to observations and is, therefore, 
relatively insensitive to the cosmological model. In any case, as was pointed 
out previously, the ISW effect is not overly sensitive to the space weighting 
function. 



2.3 X-Ray Correlations 



Given a model for the space weighting function, we then must calculate the 
fluctuations in intensity. Let A^(x, z) be the number of AGN in a volume dV 
at a position x and redshift z. If number density fluctuations are assumed to 
be related to mass density fluctuations by a redshift dependent bias b^, then 

^^(x, z) Sn{x,z) _ ^ Sp{x,z) , . s s 

= = ^^^'^^^ ^ Uz)S{^,z), (18) 

where N{z) is the average number of AGN, and 

N{^,z) = N{z) (1+ b,{z)S{^,z)). (19) 



For simplicity, throughout we will assume that the bias is independent of 
redshift. 
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The evolution of the mass density fiuctuation is characterized by, 5{x, z) — 
5(x, 0)D(2;). For flat, matter dominated cosmologies, perturbations grow pro- 
portionally to the scale factor, that is, D{z) = {l + z)~^. However, for alterna- 
tive cosmologies such as those with a cosmological constant, or with significant 
curvature, perturbations eventually cease to grow. 

The intensity of X-rays in a given direction is just the sum of the flux from 
each galaxy, 

L(n) = J dN(x, z)J^(z) (20) 



where J-'{z) is the flux we receive from a galaxy at redshift z and rn = x. 
Thus, the X-ray intensity fluctuations on the sky are 



^^(n) ^ JJ^{z)^bJz)Diz)5{x,0)dz 



(21) 



This allows one to solve for the X-ray correlation function. By Fourier expand- 
ing the present density fluctuation. 



5(x,0) = ^4e 



ikx 

k 



= Mtt E i'3i{kr)Yirr,{n^)Yi:^{n^). (22) 

k m,l 

Inserting this into equation (21), one can show that 



^ _ Am^ J J^{z)^b,{z)D{z) Ek 6^YC^{n^)jiikr 

a 



^Am'J2W)yLi^^) (23) 



where, 

,,,, = J^wf' _ 

JJ^dN 



, ^ J J^(z)^bJz)D(z)ji(kr)dz , , 



The argument deflnes a weighting function which determines how deeply in 
redshift the X-rays probe. As above, using the expression for the moments and 
integrating over the directions of k, the X-ray angular power spectrum is 

Cf = -\ k'dkP,f^{k). (25) 
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If the fluctuations have a Harrison- Zeldovich- Peebles spectrum, Pk oc k, then 
this imphes that is constant for low multipoles. The X-ray autocorrelation 
function is given by (see Eq. 2) 

Cx{e) = j^j:{2l + l)C^Pi{cos9). (26) 

Note that Cx{0) is proportional to b^. This treatment is similar to previous 
calculations of the X-ray background fluctuations (Carrera, Fabian & Barcons, 
1996; Lahav, Piran & Treyer, 1997). 

In addition to fluctuations from structure, Poisson fluctuations due to the dis- 
creteness of the X-ray sources will also contribute to the X-ray auto-correlation 
function. In general, the random placement of nearby objects causes this term 
to dominate the correlation function and the contribution formally diverges if 
arbitrarily high flux objects are considered. In practice, a high end flux cutoff 
ehminates the divergence. The Poisson contribution to the correlation can be 
written as 

{6i{^)6i{0))po^sson = J d^x'n{z')W {h, h')W {0, h')J^\z') (27) 

where W{n, n') represents the angular beam profile of the detector. This con- 
tribution to the ACF, which is independent of bias, must be corrected for in 
order to compare observations with cosmological models. 

2.4 Cosmological Models and the Power Spectrum 

The models we consider are spatially flat and are primarily composed of cold 

dark matter (CDM) with a cosmological constant. We assume the baryon 
density is constrained by Qbh'^ = 0.0125, though the dependence on baryon 
density is relatively weak. The initial fluctuations are assumed to be adiabatic 
fluctuations with a scale invariant Harrison-Zeldovich-Peebles spectrum. The 
present power spectrum is determined from the transfer function and to sim- 
plify the analysis, we have focused on models which have transfer functions 
with the same shape, as parameterized by F = Q„ihc~^^''~^^''/^"^^ — 0.25. (See, 
for example, Efstathiou (1996) and references therein.) 

We normalize the models using the variance in the CMB on 10°, as measured 
by COBE DMR (Banday, et al., 1997). From the growth of perturbations, the 
power spectrum scales approximately as (Efstathiou, Bond & White, 

1992). However, at large fl\, this dependence is softened due to the late time 
ISW effect. While a more complete analysis of the model normahzation can be 
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made using the full COBE data (Sugiyama, 1995; Gorski et al., 1995; White & 
Bunn, 1995; Bunn & White, 1996), using the 10° variance is adequate for the 
present treatment. This normalization differs from previous calculations of the 
X-ray fluctuations which normalized with respect to Ug, the standard deviation 
of the mass in a sphere of radius 8 /i^iMpc(Carrera, Fabian & Barcons, 1996). 
Roughly, the COBE normalization corresponds to a normalization of as = 
0.7 QzP'^ for the F = 0.25 models we consider here. 




k (h Mpc ) 



Fig. 2. The top figure shows the power spectrum for density fluctuations for two 
different values of F. The figure below shows where the zero lag X-ray correlation 
arises as a function of wavenumber for different resolutions. The X-ray variance, 
(dX/X)"^ , is proportional to the area underneath the curves. Dashed lines represent 
an h = 0.32, 0,^ = 1 model, while the solid lines h = 0.56, Qm = 0.5 model. 



Since we are interested in the X-ray bias, it is important to understand what 
part of the power spectrum is probed by the fluctuations in the X-ray back- 
ground. A rough estimate can be made from the flux cut of the survey and 
its angular resolution. For the given flux cut, 1 x lO^^^ergs cm~^ s~^, an AGN 
with mean luminousity L ~ 7 x lO'^'^/i^ergs s~^ would be removed from con- 
sideration if it were closer than 120/i~^ Mpc, corresponding to a redshift of 
z = 0.04. Obviously, objects dimmer or brighter than the mean luminousity 
could be closer or further away, but this distance should be a reasonable ap- 
proximation. Given a resolution of approximately 3°, this implies that the 
correlation is sensitive to structures on scales greater than 6 h~^Mpc or so. 

For a given theoretical model, this can be made more rigorous by consider- 
ing which wavenumbers contribute to the correlation. Assuming a Gaussian 



13 



window function, the smoothed correlation function can be written as, 

Cx{9) - y2{2l + l)CfPi{cos9)e-'^'+^^''Bl^'-^\ (28) 



where 9b denotes the FWHM (full width half max) resolution. Prom the above 
expression for C/, we can define 

W\9b, k)^^^ ^(2/ + l)P,ff{k)e-'^'+W^^^\ (29) 



(This is analogous to Wl^ defined by Bond & Efstathiou (1987).) With this, 
k^W'^{9B-, k) then represents the contribution to the X-ray variance for a given 
resolution from each log interval of wavenumber. 

Figure 2 shows the contributions versus wavenumber for two sample models. 
One can see that the correlations arise due to modes between 0.01— 0.5hMpc~^, 
near the peak of the power spectrum. Thus it spans the gap between those 
sampled in large scale structure surveys and those probed by CMB measure- 
ments. There are two primary effects which determine the inferred X-ray bias: 
the normalization of the spectrum and the evolution of perturbations, and 
these effects partially offset each other. The X-ray bias scales as and cr^, 
the area under the curve, scales as 



m 



2. 5 Predictions for the Cross Correlation 



Only the anisotropics created recently are correlated with the X-ray fluctu- 
ations described above. Since that is the case, anisotropics produced at last 
scattering can be ignored and only the ISW contribution to the cross corre- 
lation need be considered. The cross correlation between the X-ray and the 
microwave backgrounds is defined as 



C^^(^): 



= ^E(2^+l)^r^Kcos^). (30) 

Following the same procedures used above, these coefficients can be shown to 
be 



_i X ISW*\ 
— \"-am"-lm I 
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= edkfi{k)Pl^^ J 2^k,r)ji{kr)dr. 



(31) 



Because of the statistical nature of the background anisotropics, the corre- 
lations expressed in Equation 31 represent ensemble averages. A particular 
realization, Creaii^), of C{9) can be easily evaluted. Assuming each aim is 
gaussian distributed with variance Ci, then 



where m are numbers randomly chosen from independent gaussian distribu- 
tions with a = 1. A particular realization is, then. 



c^eim = ^ E(2^ + mrpiioos e) (33) 

where = Y.maJmO'hnI i'^^ + !)■ The distributions of C^J,(0) for 1000 re- 
alizations of several different cosmological constant models are represented in 
Figure 3. 

The spreads of the distributions reflect chance ahgnments of regions of en- 
hanced (depleted) X-ray intensity with regions of enhanced (depleted) CMB 
intensity. The latter are due primarily to the CMB anisotropics produced at 
last scattering and are, therefore, uncorrelated with X-ray emission. These 
accidental correlations can be thought of as a form of "cosmic variance". The 
curves in Fig. 3 were computed assuming full sky maps. Limited sky coverage 
increases the spread. It should be noted that widths of the distributions are 
approximately independent of A. The effect of the usual cosmic variance, i.e. 
the variance of the actual ISW effect due to the statistical nature of density 
fluctuations, is small. The proflles are roughly Gaussian and we treat this 
effect as any other source of noise in the measurement. 

Note that we have chosen to plot the quantity C-^^{9)/ax<yT which is inde- 
pendent of the bias in the theoretical calculations. All of the quantities in 
this expression are measureable. cr|. = C^(0) has been measured by COBE 
and a hmit on C^^(^) is reported in this paper. However, as discussed above, 
the measurements of a\ = C'^(O) are still ambiguous. In order to compare the 
model predictions in Fig. 3 with observations, we will take C"^(0) to be that 
predicted by a given cosmologial model modulo 6^. The results will then be 
given as a limit on the cosmological constant in terms of the unknown X-ray 
bias. 
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Fig. 3. Distribution of {XT{0)) for 1000 realizations with different values of f^A- 

3 The HEA01-A2 and COBE Maps 

The HEAOl A2 experiment was designed to measure surface brightness in 
the 0.1 — 60 keV X-ray band (Boldt, 1987). The present analysis is from 
two medium energy detectors (MED) with different fields of view (3° x 1.5° 
and 3° x 3° FWHM) and two high energy detectors (HED3) with the same 
fields of view. Counts from these four detectors were combined and binned 
in 24,576 1.3° x 1.3° pixels in an equatorial quadrilateralized spherical cube 
projection on the sky (White & Stemwedel, 1992). The combined map has an 
effective angular resolution diameter (FWHM) of 3.3° and a spectral resolution 
of roughly 2 — 10 keV (Jahoda & Mushotzky, 1989). All data used in this 
analysis were collected during the 6-month period beginning on day 322 of 
1977. 

The dominant feature in the HEAO map is the Galaxy, so all data within 20° of 
the Galactic plane and within 30° of the Galactic center were cut from the map. 
In addition, 10° diameter regions around 90 discrete X-ray sources with 2 — 
10 keV fluxes larger than 3 x 10~^^ergs s^^ cm~^ (Piccinotti et al., 1982) were 
removed. Finally, the X-ray map itself was searched for weak point "sources" 
that exceeded the nearby background by a specified amount and 7° diameter 
regions around these were removed. Cuts were made at several levels from 

4 to 10 times the photon counting noise resulting in "cleaned" maps with 
sky coverage of from 26% to 47% . The analysis that follows is for the 6 a 
cut which removed sources with fluxes greater than 1 x 10~^^ergs s~^ cm~^. 
The final "cleaned" map corresponds to 1/3 of full sky coverage; however, the 
results of this paper are largely independent of the level of the cuts. 
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Even after cleaning, the HEAO map has several components of large-scale sys- 
tematic structure which can be corrected for. If the dipole moment of the CMB 
is a kinematic effect, as it has been widely interpreted (Bennett, et al., 1996a), 
then the hard X-ray background should possess a similar dipole structure 
(Compton- Getting effect) with an amplitude of 5i/i — 4.3 x 10""^. The cleaned 
map was corrected for this effect. In addition, a linear time drift in detector 
sensivity (Jahoda, 1993) results in large-scale structure of known form. Finally, 
the 2-10 keV background shows evidence of high latitude Galactic emission as 
well as emission associated with the Superglactic plane (Jahoda, 1993). We 
modeled these latter three effects as a combination of a hnear time drift, a 
Galactic secant law, the Haslam 408 MHz Galactic radio map, and a simple 
pancake Supercluster. This model was linearly regressed to the cleaned data 
and subsquently subtracted from the map. Of the four parameters, the time 
drift and secant law were most significant. Correcting for these effects signifi- 
cantly reduced the large scale structure in the X-ray autocorrelation function 
but had little effect at small angular scales as expected (see Section 3.1). When 
the Compton-Getting dipole was included in the fit (3 additional parameters), 
the results did not change significantly. 

Because of the echptic longitude scan pattern of the HEAO sateUite, sky 
coverage and, therefore, photon shot noise were not uniform. However, the 
mean variance of the cleaned, corrected map, 1.63 x 10^^ TOT counts/sec, 
is considerably larger than the mean variance of photon shot noise, 0.67 x 
10~^ TOT counts/sec, where TOT counts/sec is the standard HEAOl A2 
normalization for the 2-10 keV band (Allen, Jahoda & Whitlock, 1994). The 
correlation analysis of Section 3.1 demonstrates that the variance in the HEAO 
map is due to small angular scale (< 5°) intensity fluctuations; therefore, 
the X-ray map is dominated by "real" (not photon noise) structure. For this 
reason, in the correlation analyses that follow, we chose to weight each pixel 
equally. 

The CMB map was constructed from the 53 GHz and 90 GHz 4-year COBE 
DMR maps as obtained from the National Space Science Data Center (Bennett, 
et al., 1996b). Each map consists of 6144 2.6° x 2.6° pixels in an ecliptic quadri- 
lateralized spherical cube projection, i.e. half the resolution of the X-ray map. 
The 31 GHz maps have considerably larger instrument noise and Galactic 
contamination and were not used in this analysis (Bennett, et al., 1996a). 
These four temperature maps (A and B channels for each frequency) were 
converted from antenna to thermodynamic temperature and then combined 
in a straight average to form the composite CMB map. This straight average 
map has somewhat larger instrumental noise than a noise weighted average 
map; however, the noise in the crosscorrelation function is dominated by "cos- 
mic variance" (see Section 3.1) and we felt that a straight average would be 
more likely to minimized unknown systematic effects in the composite map. 
The same Galatic cut (within 20° of the Galactic plane and within 30° of the 
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Galactic center) was applied to the CMB map as to the X-ray map which 
results in 64% sky coverage. 

The only correction made to the cleaned composite CMB map was to fit and 
remove the dipole moment. Since the quadrupole moment was included in the 
analysis of section 11-C, no attempt was made to remove it from the data. 
The removal of a secant law fit to high Galactic latitude emission made no 
significant difference in the cross-correlation function (see Section 3.1) and so 
the GMB map was not corrected for Galaxy emission. In any case, the analysis 
of this paper concerns the shape of the cross-correlation function at relatively 
small angular scales (< 20°) and is not overly sensitive to the presence of large 
scale structure in the map. 

As with HEAO, the GOBE sateUite sky coverage was not uniform; therefore, 
instrument noise per pixel is also not uniform over the sky. In this case the 
mean variance in the cleaned, corrected GMB map is 5.0 x lO^^mK^ whereas 
the mean variance of instrument noise is 4.1 x 10~^mK^, i.e. at the resolution of 
the map, instrument noise dominates real structure. However, when smoothed 
in 10° bins, instrument noise variance becomes 3.4 x lO^^mK^ compared to 
the 9.0 X lO-'^mK^ fluctuations of the GMB (Banday et al., 1996). Therefore, 
in the correlation analyses that follow, we also chose to weight each GMB pixel 
equally. 

3.1 The X-ray Auto- Correlation Function 

The intensity auto-correlation (AGF) is calculated using, 

{XX{d))=Y:{i,-l){i,-l)/Ne, (34) 

where the sum is over pair of pixels separated by 9 and is the number of 
such pairs. Note that the pairs are given equal weight as discussed above. 

The results for the cleaned X-ray map are shown in Fig. 4. The value at ^ = 
has been corrected for contribution of uncorrelated photon shot noise. To check 
for the contribution to the AGF of weak X-ray sources in nearby galaxies, 4° 
diameter holes around all (2367) galaxies in the Tully Atlas of Nearby Galaxies 
(TuUy, 1988) were cut from the map. These cuts did not significantly change 
the AGF. 

Taking into account the finite angular resolution of the HEAO map, the ampli- 
tude of the AGF is somewhat smaller than that found by Soltan et al. (1996) 
for the ROSAT 0.5 — 2 keV. However, it should be noted that the Soltan et al. 
results are inconsistent with other measurements (Ghen et al., 1994; Garrera 
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Fig. 4. The auto-correlation function for the HEAOl-Al X-ray map, normahzed in 
terms of the mean intensity. Shown for comparison is the theoretical prediction for 
a COBE normalized CDM model with an effective smoothing of 3.6° FWHM. 

et al., 1993) leading them to hypothesize diffuse emission from large (20 Mpc) 
clouds of gas. 

It should be emphasized that the ACF in Fig. 4 has not been corrected for 
Poisson noise due to unresolved sources. The model of Comastri et al. (1995) 
indicates that source confusion is substantially less that the measured value 
while the source counts of Piccinotti et al. (1982) and Butcher (see (Carrera 
et al., 1993)) imply a substantial amount of source confusion. On the other 
hand, poor angular resolution results in shape of the ACF being dominated 
by the beam profile. For these reasons we have not chosen to estimate the 
correction for source confusion and the data in Fig. 4 should be considered an 
upper limit to the ACF due to clustering of X-ray sources. 

Fig. 4 also shows the theoretical ACF arising from a CDM power spectrum 
with A = 0.5, (r = 0.25), a bias of bx = 3.9, and smoothed by a 3.6° gaussian 
beam. Recall that the height of the theoretical ACF is proportional to 6^. The 
sources were assumed to be distributed as in the model of Comastri et al. 
(1995) with a flux cutoff equal to that applied to the real data. Clearly this 
bias is a good fit to the model; however, to the extent that the observed ACF 
is contaminated by source confusion, it should be considered an upper limit 
for this model. Raising F has the effect of adding small scale power, and so 
would reduce the bias required to fit the X-ray ACF. 
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3.2 The HEAO/COBE Cross Correlation 



The cleaned, corrected X-ray map was reprojected onto the sky to match the 
COBE map, i.e., an echptic quadrilaterahzed spherical cube projection with 
2.6° X 2.6° pixels. Let Wj be the number of original pixels combined to form one 
pixel in the new projection. The value of Wj for most pixels was 4; however, 
because of the transformation to ecliptic coordinates, Wj ranged from 1 to 6. 
The terms in the cross-correlation function (CCF) are weighted with these 
values. The cross correlation is then computed as. 



{Tx{e)) = J2 x.^iTj/ E 



(35) 



where the sum is over all pairs with angular separation 0. 
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Fig. 5. The cross correlation between the four year COBE map and the HEA01-A2 
X-ray map, normalized by the rms fluctuation of each map. Also shown are the 
expected correlations for a range of values for the cosmological constant. (The stan- 
dard CDM model predicts no correlation.) The bold lines show the rms correlations 
resulting from noise and cosmic variance. 

The resulting cross correlation is shown in Fig. 5, where the normalization is 
with repect to the measured rms fluctuations of the two maps. Also shown 
are the theoretical predictions for cosmological constant models with varying 
A. As described previously, the Hubble constant has also been varied to keep 
the shape parameter, F, flxed. Finally, the solid curves represent the la un- 
certainty due to the combination of cosmic variance, instrument noise, and 
photon shot noise. The spread in the distribution illustrated in Fig. 3 did not 
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include the effects of map cleaning, map correction, and reduced sky cover- 
age. In addition, there are quite likely additional fluctuations in the X-ray map 
that do not appear the X-ray model, e.g. relatively nearby clusters of galaxies. 
Therefore, we chose to compute the total "noise" by a Monte Carlo method 
wherein the real, cleaned X-ray map was cross-correlated with an ensemble of 
400 randomly generated CMB maps. These latter maps consisted of random 
realizations of instrument noise and cosmic structure which was constructed 
so as to have the same ACF profile as the actual COBE data. All the Monte 
Carlo CMB maps were cleaned (Galaxy cut) and corrected (dipole removal) 
with the same procedure as used on the real map. The resultant total noise is 
correlated on angular scales of about 10°. This is due primarily to the finite 
angular resolution of the COBE map. It is clear from Fig. 5 that the data are 
consistent with the hypothesis that the two maps are uncorrelated. 

The above results assume that fiuctuations in the X-ray map are dominated 
by Poisson fluctuations of unresolved sources. If there is a significant contribu- 
tion to the X-rays from this Poisson noise, then both the inferred X-ray bias 
and the predictions for the theoretical correlation are reduced by a factor of 
(7 structure/ {cT Structure + (^Poisson)^ ■ Therefore, for any particular bias, the model 
amphtudes should be scaled as bx/5.6 O,'^^. In any case, it seems unhkely that 
the bias would be less than unity. 

The upper limit of the CCF at ^ = is consistent with previous correlation 
analyses (Boughn & Jahoda, 1993; Banday, et al., 1996) between HEAOl- 
A2 and temperature anisotropy maps. The correlation analysis of the COBE 
DMR maps and the ROSAT PSPC All-Sky Survey (Kneissl et al., 1997) shows 
somewhat more structure than Fig. 4. This is presumably due to a larger 
Galactic contriubtion at 1 keV. 



4 Discussion 

As previously discussed, the observed HEAO/COBE cross-correlation is con- 
sistant with there being no underlying ISW effect, i.e. no cosmological con- 
stant. Since the effective errors computed by Monte Carlo analysis (see Section 
3.2) are consistent with being derived from a correlated gaussian distribution, 
we can quantify this by performing a maximum likelihood analysis. The like- 
hhood is defined by 

C = (27r)-^/2[det R,j]"^ exp[-^i?-.^(y, - y,){yj - y,)], (36) 

where Rij is the covariance matrix and tji is the theoretical prediction for the i**^ 
observation. To a good approximation, the theoretical cross-correlation curves 
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in Fig. 5 have the same shape and only differ in amphtude (which depends 
only on Qa), 



{TX{0))^A{A)m 



(37) 



where /(0°) = 1. Maximizing the likelihood for this amplitude, one finds, 



4 = Y.^ij'R'uiRjififk + 3 permutations )/(2 ^ i^r.^/./.j^. (39) 



Since "cosmic variance" is the dominant uncertainty, we calculate the covari- 
ance matrix from the correlations of the X-ray map with the simulated CMB 
maps. Although this implicitly assumes that there are no correlations between 
the maps, it is a good approximation when the ISW efi^ect is small as discussed 
in Section 2.5. 

The most likely amplitude is found to be A = —0.008 ± 0.042, based on the 
measured correlation for 9 < 25°. (At larger separations the models predict 
little correlation.) This corresponds to a 95 % CL upper limit oi A < 0.061 
(1.65(7), and a 98 % CL upper limit oi A < 0.078 (2.05a). The amplitudes. A, 
for different models are summarized in Table 1 along with the X-ray biases 
inferred from the X-ray map assuming all of the X-ray fluctuations result from 
large scale structure. Under this assumption, > 0.5 is just ruled out at the 
95 % CL. As discussed above, if there is a significant contribution to the X- 
rays from Poisson noise, then both the inferred X-ray bias and the predictions 
for the theoretical correlation are reduced by a factor of cr structure/ {^structure + 
Poisson)'^ ■ Included in the table are the expected amplitudes for a fixed X-ray 
bias, br^ = 1 or bx = 2. If b^ = 2, Qj>^ > 0.7 is still ruled out at the 95 % CL. In 
the most pessimistic case, with b^ = I, none of the models can be ruled out. 

The analysis we have presented has been conservative in the sense that we have 
only used the X-ray temperature cross correlation to constrain the models. For 
nonzero cosmological constant, the spectrum Ci of temperature anisotropics 
shows a rise at low / which is absent in the COBE spectrum. Inclusion of 
the latter in our analysis would therefore strengthen our bound on (For 
example, see Sugiyama & Silk (1994).) We have chosen not to do so because 
the low / rise is a more model-dependent effect than the X-ray temperature 
correlation, it being possible to remove the former completely by a tilt in the 
primordial power spectrum. Nevertheless, a future analysis of a wider class of 
models should include the Q information as well. 



(38) 




and the variance in this amplitude is given by. 
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Table 1 

The inferred biases and the expected cross correlation amplitutudes, 
{XT{0°))/axO'T for a range of models with F = 0.25. The inferred biases 
and the first column of amplitudes assume that the X-ray fluctuations result solely 
from large scale structure. The other amplitudes assume a fixed X-ray bias, bx = I 
or bx = 2. 

It is important to search for as many independent constraints on a cosmologi- 
cal constant as possible, because each is subject to different kinds of systematic 
errors. At present, the best established limits on A come from gravitational 
lensing studies (Maoz & Rix, 1993; Kochanek, 1995). Kochanek (1995) finds 
that the lack of observed lensing events implies Q\ < 0.66 at the 95% confi- 
dence level. In addition, recent studies of the deceleration parameter measured 
in supernovae searches have placed an upper limit of Q\ < 0.51 at the 2a level 
(Perlmutter, et al., 1996). 

The limit on the cross-correlation amplitude can also be used to constrain open 
models, where the correlations could be much higher for a given Q^, if the 
X-rays probe at high enough redshifts (Kamionkowski, 1996). The constraints 
on open models thus are likely to be much more sensitive to assumptions 
of the luminousity evolution of the sources than the cosmological constant 
constraints. 



5 Conclusions 

Our primary aim in this paper has been to demonstrate how correlations, 
or the absence of correlations, between the microwave background and deep 
probes of structure can lead to constraints on cosmological models. The sur- 
veys we used were not ideal for this purpose, but even so they resulted in a 
relevant constraint on the cosmological constant. 

Surveys with better angular resolution will improve the signal to to noise 
and thus the confidence with which we can rule out a particular cosmology. 
With degree scale resolution in both the temperature and X-ray maps, the 
signal to noise could be improved by 50%. For the hnear effect which we are 
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considering here, little is gained by looking at maps with resohition on scales 
smaller than a degree or two. The CMB photons do not have time to receive a 
substantial shift from the time varying potential when passing through smaller 
fluctuations. Since the correlation is restricted to large scales, full sky maps 
are essential to reduce uncertainties due to accidental correlations. 

Future satellite missions, most notably MAP and the Planck Surveyor, will 
provide full sky maps with sufficient resolution (0.3° and 0.1°, respectively) to 
be explore this effect fully in the CMB sky. Data from these probes should be 
available in the next 5-10 years. The prospects for better data for the X-ray sky 
in the near future are uncertain. The ROSAT survey is sensitive only at lower 
energies (< 2 keV) and so is likely to be contaminated by galactic and other 
foregrounds. To search for this effect, ideally one needs a full sky survey of the 
hard XRB (2-10 keV) which has degree scale resolution and which is sensitive 
enough to allow an unambiguous subtraction of Poisson fluctuations. Such a 
survey would also provide a measurement of structure on scales larger than 
is possible in current optical surveys and would fill in the gap in the power 
spectrum which lies between the galactic surveys and the measurements of 
the microwave background. This kind of mission is now being considered (K. 
Jahoda, private communication). 

One can also consider looking at other possible probes of structure at high red 
shifts. The formalism for calculating the correlations presented here for X-rays 
transfers easily to other kinds of probes. Possibilities include surveys of radio 
galaxies or quasars, or other objects which probe to large redshifts. Galaxy 
surveys such as the Sloane DSS, which should probe to 2; ~ 0.5, might also 
be sensitive to this effect and have the added benefit of containing redshift 
information for the individual sources. An advantage of these probes is that 
their large scale clustering properties are (or will be) known and the bias 
implied by a given model can be easily computed. However, if the surface 
density of these objects is too low (e.g. quasars) then the statistical noise in 
the cross-correlation function will be large and the effects discussed in this 
paper will be unobservable. 

In this paper we have focused on testing for a cosmological constant, but in 
fact this is a rather indirect means of testing for A. The results we quote are 
based on a fairly specific model, one which has cold dark matter as well as 
adiabatic, scale invariant initial fiuctuations. It is important to realize however 
that rather than being sensitive merely to one specific cosmological model, 
the effect we have considered occurs in any cosmology where the gravitational 
potential has evolved recently. This includes models that are open or where 
structures arise from cosmological defects and we are presently working to 
understand this constraint in these other contexts. 

Whatever the cosmological model, the ISW effect reflects important funda- 
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mental information relevant to structure formation. The presence of large, 
empty voids, the apparent lack of mergers and the shapes of clusters of galax- 
ies provide evidence, tentative as yet, that the process of structure formation 
might be slowing down. If this is the case, then the decay of the linear gravita- 
tional potential makes a late time ISW effect inevitable and offers an important 
observation to confront these other clues about how structures formed in the 
universe. 
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